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An exact expression for the cross polarization between two spin- 1/2 particles is derived from the 
quantum Liouville equation. This is given in the form of two integrodifferential equations. These 
can be solved exactly in the static case (no sample spinning) and a powder average easily performed 
numerically. With magic angle spinning, the neglect of certain interference terms simplifies the 
numerical calculation. A further assumption decoupling the calculation of the sidebands gives a very 
simple formula that is capable of giving a qualitative interpretation of all experimental observations. 
Examples are given illustrating typical buildup curves and CPMAS matching profiles. 
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I. INTRODUCTION 



The most frequently used technique to improve the 
sensitivity of nuclei with low gyromagnetic ratios and/or 
low isotopic abundances in solid-state NMR is the 
Hartnmnn-Hahn cross polarization (CP) in the rotating 
frameEJ. With this technique, polarization can be trans- 
ferred from sensitive and often abundant nuclear spins 
/ (e.g., protons) to the insensitive nuclear spins S (e.g., 
13 C, 15 N). It is driven by the flip-flop transitions where 
both dipolar coupled spins / and S change their mag- 
netic quantum number in an energy conserving process. 
Since this transfer of polarization takes time to develop, 
CP can be observed only if the whole spin system has 
reasonably long relaxation times. The CP spin dynamics 
has already been discussed more then 25 years ago and 
since then has been the subject of a variety of additional 
studies. The polarization of the abundant I spins (usu- 
ally with a high gyromagnetic ratio) is spin-locked in the 
rotating frame by a radiofrequency (rf) field. When a 
second rf field of appropriate magnitude is then simulta- 
neously applied to the dipolar coupled rare S spins (usu- 
ally with a low gyromagnetic ratio), flip-flop processes 
become energy conserving and thus allow the transfer of 
polarization from the / spins to the unpolarized S spins. 
An overview of the Cross Polarization under Magic Angle 
Spinning (CPMAS) can be found in a number of places, 
a partial list of which is Refs. ||^,||,[5]. The emphasis 
in this paper is on an exact theoretical formulation of 
CPMAS and no attempt is made to review the experi- 
mental methods and/or results that have been achieved 
using them. 

Several theoretical approaches have been presented for 
the description of CPMAS. One such approach is based 
on there being a spin temperature, which assumes the 
spins are in thermal equilibrium in an appropriate rotat- 
ing frame, the recent paper by MeierQ exemplifies how 
this approach can be used. Another method, due to 
Miillcr et alEl, resulted in an empirical fitting formula 
which has been extensively used by many workers. A 
third teflpe of approach is based on the use of floquet 
theoryBIa and_a comprehensive program GAMMA has 
been writtenE3 which can be used for treating spin dy- 
namics and numerically solving the quantum Liouville 



equation. A method of describing the spin dynamics 
in terms of rotations, by parts of the_hamiltonian has 
been used by LevittU and Zhang et altL3. Finally, there 
are methods that use the Moril3 approach to nonequilib- 
rium statistical mechanics. One method is based on pro- 
jecting onto selected parts of the hamiltonianEJ whereas 
anothenI3 uses standard approximation methods for esti- 
mating the memory kernel. In contrast to this variety of 
methods, the approach presented here utilizes the special 
nature of the hamiltonian in the doubly rotating frame, 
and the specific properties of spin- 1/2 systems to obtain 
an exact integrodifferential equation (correctly a pair of 
such equations) for the cross polarization. A solution of 
these equations is complicated by any time dependence 
of the dipolar coupling, while an exact solution is imme- 
diately available for a time independent and orientation 
fixed coupling constant. The powder average of the re- 
sulting solution can be done numerically. 

While the spin temperature approach is very appeal- 
ing in that it provides a relatively simple explanation for 
the qualitative behaviour of cross polarization, it has a 
number of weaknesses. First of all, it is noticed that the 
spin system suffers no (or at most negligible) relaxation 
during a cross polarization experiment. Thus there is 
no basis for an equilibrium to be obtained and thus no 
basis for assuming a thermal distribution. An argument 
might be made for ergodic behaviour but we do not find 
that appealing. Second, the thermal equilibrium is to be 
in the doubly rotating frame, in which the hamiltonian 
must still be time dependent if magic angle spinning is 
performed. Thus a selection has to be made in some way 
to extract a time independent hamiltonian in order to 
have a thermal distribution. Exactly how this is done is 
arbitrary. Thirdly, only one Hartmann-Hahn matching 
condition can be described at a time, after selecting out 
the right time independent hamiltonian. The theory pre- 
sented here is a purely dynamical theory (no relaxation 
is involved), is capable of treating the magic angle spin- 
ning exactly (approximations simplifying the calculation 
are made, but all time dependence could be included), 
and all four peaks of the CPMAS matching profile are 
obtained from the same calculation. 

The paper starts by reviewing some relevant properties 
of spin-1/2 systems together with their rotational prop- 
erties. These are used in Sec. Ill to examine the specific 
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nature of the spin hamiltonian in the doubly rotating 
frame and thus show that the relevant operator space 
for describing cross polarization is only 6 dimensional. 
Section IV then derives a formal expression for cross po- 
larization, described in terms of two integrodifferential 
equations. Section V gives an exact expression for the 
static (no spinning) case and exemplifies the purely os- 
cillatory time dependence for a constant dipolar coupling 
strength (fixed direction of the vector joining the parti- 
cles containing the spins), and also the destructive inter- 
ference effects that occur in a powder sample. Sec. VI 
describes the effect of magic angle spinning and a short 
discussion ends the paper. Figures illustrating a variety 
of different behaviors are presented. 

II. A SYSTEM OF TWO SPIN-1/2 PARTICLES 

The spin operators for the two particles are labelled as 
S x , S y , S z and I x , I y , I z . There are also the identities 
I5 and 1/ for the respective spin spaces. As spin-1/2 
particles, these operators satisfy the product conditions 

S x S y = 2*-' z > Ixly — 2^ z ' W 

S X S X — -Is, I x Ix — ^1/, (2) 

and their cyclic permutations. The first set of product 
conditions are of course just part of the commutation 
relations 

[S x ,S y ] — = iS z \I x ,Iy\— — Hz, (3) 

and their cyclic permutations. [Throughout this work, 
h is set equal to 1.] The four operators for each spin 
(including the identity) implies a total operator space 
of 16 dimensions, which exactly coincides with the di- 
mension of a matrix representation of a system of two 
spin-1/2 particles. Thus any operator can be expressed 
in terms of the sum of products of S and / spin operators. 
For simplicity of presentation, the identity operators are 
usually not written explicitly. Thus the operator S z is 
understood to stand for the product S z li, etc. 

Another important property of spin operators (not just 
of spin-1/2) is that the commutator of a spin operator 
generates a rotation. Setting the commutator superop- 
erator 

S Z A={S Z ,A}_, (4) 

the exponential superoperator 

e - ws * S x = e~ ws * S x e ies ° = S x cos 9 + S y sin 9 (5) 

rotates the S'-spin operator on which it acts (here S x ) 
about the z axis. X z is used to denote the corresponding 
commutator superator for l z . Analogous commutator su- 
peroperators for the other spin directions can be defined 
but are not needed in this work. 



III. NATURAL PARTITIONING OF THE 
OPERATOR SPACE IN THE 
DOUBLY-ROTATING TILTED FRAME 

This is phrased in the tilted frame, namely the rf fields 
are labelled as being in the z direction and the spin op- 
erators of the truncated (often referred to as the secular 
part of the) dipolar hamiltonian lying along the x direc- 
tion. Thus the hamiltonian isE3 

H(t) = u u I, + lu 1s S z + 2b(t)I x S x . (6) 

The angular velocities u>u = \jiBu\ and luis = \jsBis\ 
are proportional to the radio frequency field amplitudes 
Bu and Bis- In the laboratory frame, these rf fields 
are oscillating at the respective resonance Zeeman fre- 
quencies determined by a large static magnetic field in 
the x direction. b(t) is the (possibly) time dependent 
heteronuclear dipolar coupling strength. Fundamentally 
this coupling strength depends on the angle 9 between 
the vector r pointing from the nucleus having its spin 
labelled / to the nucleus whose spin is labelled S relative 
to the Zeeman field direction x, being explicitly given by 

b{6) = b Q (cos 2 9~ 1/3), (7) 

where bo = — 3/io7/7s , ?i/(87rr 3 ) is the dipolar coupling 
constant. Ignoring any internal (rotational or vibra- 
tional) motion of the molecules in the chemical system, 
a time dependence arises only through the physical spin- 
ning of the sample about an axis at the magic angle 
[arccos(l/-\/3)] to the Zeeman field (a; direction). Thus 9 
becomes time dependent and b(t) is written as a short- 
hand for b(9(t)). There is inherently a remaining orien- 
tation dependence of b(t), associated with the original 
orientation 9(0) of r. In a powder sample there is a ran- 
dom distribution of such angles so an average over 9(0) 
needs to be made for such samples. But for the analysis 
of the spin dynamics, only the possible time dependence 
of b(t) is important, so only this time dependence is ex- 
plicitly displayed. Later on, when a powder sample is 
considered, the extra averaging over initial orientations 
(or its equivalent) is made. 

It is also noted that, if b is time independent (and has 
a fixed interspin orientation, 9), the energy eigenvalues 
of this hamiltonian are ±(l/2)^/6 2 + (cjj ± lus) 2 - These 
eigenvalues have been noted (at least) by Levitt, et allij, 
but their method of solving the spift-jdypamics, as do 
some other treatments, for exampletJIla'tZl, is based on a 
separation of the hamiltonian, Eq. (||), into two commut- 
ing parts. Neither the eigenvalues, nor their respective 
eigenvectors, seem to be standardly quoted in the liter- 
ature. They are listed for reference in Appendix A. The 
method used in this paper does not explicitly use these 
eigenvalue properties nor the above noted separation of 
the hamiltonian into commuting parts. 

The purpose of this section is to catalog how the 16 
spin operators behave under the action of the commuta- 
tor superoperator C(t) for the hamiltonian, 

C(t)A = [H(t), A]_ = C I{ A + C d (t)A. (8) 
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This is listed as a sum of radio frequency and dipolar 
parts. It is to be shown that this set of operators is par- 
titioned into four dynamically separate sets of operators. 
Crucial to this partitioning is the vanishing of commuta- 
tors of the form 

\S X I X i Syly] — = S x [4 , Syly] — ~\~ [S X: Syly] — I x 
= S X Sy(H Z ) + iS Z IyI X 

= (i/2)S z (U z ) + iS z {-i/2)I z = 0. (9) 

This relation is special to spin- 1/2 systems since it de- 
pends on the special product conditions, Eq. ([!]), of the 
spin-1/2 operators. 

Starting with the operator 4, the action of C(t) on 
successively appearing operators is 



CI Z = 



2ibS x Iy i 



£S X Iy = iUJ 1S SyIy — i^llS X I X + ~ M Z 



CSyly 

CS X I X 

£SyI X 

cs z 



Jy±y ^.. ±J ^^^^ , 

= —iuisS x Iy — itonSylx 
= iuisSylx + iuJiiS x I y 

= -iuisS x I x + iuiuSyly + -bS, 
— - *2/zbSyJ x . 



(10) 



The possible time dependence of C and b has not been 
indicated in these equations. This gives a closed set of 
6 operators, namely I Z ,S Z , S X I X , S x I y , S y I x and S y I y . It 
is this set of operators that govern the behaviour of a 
CPMAS experiment. The set of commutators ( |l(i| ) have 
some special properties which are useful for the subse- 
quent analysis. Specifically it is noticed that the dipolar 
hamiltonian acts only to couple the single 1-spin oper- 
ators S z and I z to the product spin operators, whereas 
the rf operators act only on the product spin operators. 

Another closed set consists of the identity lsl/ as its 
sole member, since this operator is unchanged by C. So 
it is a set all by its own. The same is true for the product 
operator S Z I Z since it commutes with both rf operators 
and, because of the commutator type of Eq. (^|), it also 
commutes with the truncated dipolar operator. Finally 
there are the remaining 8 operators, which can easily 
be shown to be all dynamically coupled, so there is no 
further partitioning of the spin space. 

Since the 16 operators span the operator space, any op- 
erator, and in particular the density operator p, can Jae 
expressed as a linear combination of these operatorsEa. 
For the CPMAS set of operators, the relevant space con- 
sists only of 6 operators plus the identity, the latter being 
retained in order that the density operator be normalized, 
namely that the trace over both spins (a total of four 
Dirac states) is consistent with Trp = 1. Of these opera- 
tors, only the identity has a finite trace, so p must have 
(l/4)lsl/ as one of its terms. Recognizing the expansion 
coefficients as proportional to the corresponding expec- 
tation values, the CPMAS density operator can thus be 
expressed as 

p = (1/4)151/ + S Z {S Z )+ 4(4) + AS X I X {S X I X ) 



+ AS X Iy(S X Iy) + 4SyI X (Sylx) + ASyIy{SyIy). (11) 

Thus the study of the evolution of the density operator is 
equivalent to the study of the evolution of the expectation 
values of these observables. 

It is appropriate to subdivide this set of operators fur- 
ther, namely into the set of 1-spin operators S z and 4, 
plus the identity to preserve normalization, and the re- 
maining set of 2-spin operators. Appropriate for this 
subdivision, a projection superoperator V is defined that 
selects out the 1-spin part of the operator space, thus 



V = (1/4)1 S 1/Tr + S z TrS z + 4Tr4 



(12) 



where the trace is to be taken over both spin spaces. 
Necessarily, this superoperator is idempotent TV = V, 
which determines the numerical factors. 1 — V is written 
for the projector onto the remaining (2-spin) operators. 
Because of the commutation relations, it is noticed that 



TCP 

(1 - V)C(1 - V) 
VC{1-V) 
(l-V)CV 



= o, 



£ r f — uiisS + will, 

vc d , 



(13) 



separates the action of L within the V and 1 — V spaces 
to be entirely due to the rf fields (vanishing for the V 
space), while the coupling between these two subspaces 
is due solely to the dipolar hamiltonian. The retention 
of one projector in each of the subspace coupling cases 
serves to select which space is to be acted upon and which 
is to be mapped into. Again it is emphasized that any 
time dependence of the dipolar coupling strength has no 
effect on this separation of the different parts of C. 



IV. THE SPIN OPERATOR DYNAMICS FOR 
CROSS POLARIZATION 

After a 90° pulse to align the /-spin along the z axis 
(doubly rotating tilted frame), the time evolution of the 
expectation value (S z ) of the S'-spin is monitored exper- 
imentally. Thus the calculation is to use the value of 
(4)(0) as the only nonzero spin expectation value as the 
initial condition and to calculate the consequent time de- 
pendence of (S z )(t). Since it is the expectation values 
of S z and 4 that are of importance, it is the partition 
consisting of these two operators that is relevant for the 
analysis of these experiments. The projection superop- 
erator V defined above carries out this selection. 

The density operator p for the spin system satisfies the 
quantum Liouville equation (von Neumann equation) 



(14) 



Again Ti is set to 1. The expectation value for any ob- 
servable A is calculated as the trace (A) = TrAp whose 
time dependence can be calculated according to 



d(A) 
dt 



TtA 



dp 



-iTrACp = iTr(£A)p, 



(15) 
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in the Schrodinger picture (first form) or Heisenberg pic- 
ture (last form). From the commutation relations (10), 
it is a set of 6 coupled equations that must be solved 
in order to predict the time evolution of (S z )(t). Since 
only the S z and I z part of the operator space reflects the 
experimentally observed quantities, a separation of the 
dynamics to emphasize_ihese can be useful. The gener- 
alized master equatiorO approach does this in a natural 
way. This method is summarized in the following sub- 
section. 



A. Derivation of the generalized master equation 

The derivation of the generalized master equation is 
reviewed in its generality since there is no immediate sim- 
plification of the derivation by specializing the superop- 
erators to the present case. Note that a time dependence 
of £(t) is allowed for here, which is standardly not con- 
sidered. 

Define the relevant part of the density operator as p\ = 
Vp and the irrelevant part as the remainder p2 = (1 — 
'P)p. Then these parts evolve according to 



VC(t)p(t) = VCVpi (t) + VC{t){\ - V)p 2 (t), 



.dpi(t) 

= (i-w)p(*) = (i-nc(i-n*(*) 



-(l-V)C(t)V Pl (t). 



(16) 



Here the possible time dependence of the coupling super- 
operators VC(t)(l - V) and (1 - V)£(t)V are explicitly 
displayed. Solve formally for p%(t), 



P2(t) = e 



-i{\-V)C{\-V)t 



p 2 (o) 



-lie 
lo 



(l-V)C(t')V Pl (t')dt'.(17) 



Note that it is assumed that there is no time dependence 
of (1 — V)C(1 — V), so the exponential evolution super- 
operator is correct as written. This is the case for the 
present problem. Substitute this equation for p2(t) back 
into the differential equation for pi(t), 



i% = VCp = VCV P x 
at 

+VC(l-V)(t)t 



-i{l-V)C{l-V)t 



P2(0) 



x{i-r)c{t')v Pl (t')dt'. (is) 

This is the generalized master equation of ZwanzigB, 
with an additional time dependence of the coupling su- 
peroperators. It describes the exact evolution of the rele- 
vant observables while taking into account the role of the 
irrelevant observables in this evolution by means of an 
integral over the time history (memory) of the relevant 
observables. 



B. Application of the generalized master equation 

For a cross polarization experiment involving spin- 1/2 
particles, the special properties, (|l3|), of the Liouville su- 
peroperator allows all projection superoperators to be 
dropped in the generalized master equation for pi(t). 
Moreover, the initial condition on the density operator 

p(0) = Pl (Q) = (1/4)151, + I z (I z )(0) (19) 

lies entirely in the relevant space, so the generalized mas- 
ter equation reduces to 



= -/ dWe-^-V&WpxWde. (20) 



dpijt) 
dt 



This is exact for this system! 

On taking the trace of this equation with S z and using 
the expansion 

pi(t) = Vp(t) = (1/4)131/ + S z (S z )(t) + I z (I z )(t) (21) 

for pi (i), the equation of change for the expectation value 
(S z )(t) is calculated to be 



d(S z )(t) 
dt 



b(t) [K ss (t-t')(S z )(t') 
+K SI (t-t')(I z )(t')}b(t')dt'. (22) 



Here the (possibly time dependent) dipolar coupling 
strengths are written out explicitly and it should be noted 
how they depend on different times. The memory kernel 
Kss(s) is defined as 

Kss(s) = 4:T±[S z ,S x I x }-e~ i{u)lsS * +UllX * )s [SxIx,Sz}- 
= 4Tr{i)S y I x e- l ^ sS *+^< x ^ s (-i)S y I x 
= 4TrS y I x [S y cos(loi S s) - S x sin(wi S s)] 

x [I x cos(o;i/s) + Iy sin(wi/s)] 
= cos(ujiss) cos(o;i/s) 

= [cos((cJi S + uj u )s) + cos((wi5 - u u )s)] . (23) 

In the first term the commutator [S e ,S x I x ]- arises by 
using the invariance of the trace to interchange the order 
of commutation arising from TrS z £d ■ ■ ■. The memory 
kernel Ksi(s) is calculated in a similar manner to give 

K S i(s) = (1/2) [cos((wi S + uj u )s) - cos((wi S - uxi)s)] . 

(24) 

The corresponding equation for (I z )(t) is obtained by 
interchanging S and / labels on the various terms. Specif- 
ically this gives 



d(h)(t) 
dt 



b(t) [K IS {t-t')(S z ){t') 
+K H (t-t')(I z )(t')]b(t')dt'.(25) 
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The memory kernels 



and differentiating again 



K H (s) = 1/2 [cos((wis + uj u )s) + cos((wi S 
Ki S {s) = 1/2 [cos((wi S + Wu)s) - cos((wi5 



wu)s)] , 
wu)s)] , 
(26) 



are of the same form as the previous pair. As a conse- 
quence, the equations for (S z ) and (I z ) can be simplified 
when rewritten in terms of the sum and difference expec- 
tation values, namely 



d(S,-I,)(t) 



dt 



and 



b(t)b(t')cos(A(t-t'))(S z -I z )[t')dt' 

(27) 



d(s z + i z ) (t) = r* b ^ b ^ cos( x( t _ t ^( Sz+Iz) (t / )rft / ) 

Jo 

(28) 

u>li and sum S = 



dt 



Wis 



involving the difference A 
uj\s + u>u frequencies. 

It is thus shown that the calculation of the cross polar- 
ization separates into two independent parts. This com- 
pletes the part of the calculation which involves the spin 
operators, the remaining part of the calculation requires 
solving the two integrodifferential equations. In order to 
carry this out, it is necessary to specify the detailed time 
dependence of the dipolar coupling strength b(t). The 
static [b(t) time independent] and magic angle spinning 
cases are discussed in turn. 



V. STATIC CROSS POLARIZATION 

In the static case, the dipolar coupling strength b is 
time independent, but still depends on the angle of the 
vector r separating the S and / spin particles and the 
direction (x) of the static magnetic field. In a powder 
sample there is a random arrangement of this direction so 
an ensemble average of (S z )(t) over 8 needs to be taken. 
But this is to be done after (S z ) has been calculated for 
each fixed direction. 



A. Fixed orientation of the dipolar coupling 
strength 

The expectation of the difference of the spins D(t) 



(S z — I z )(t) is considered first. According to Eq. ( J2T[ ) 
this satisfies the integrodifferential equation 



dD(t) _ h2 
dt 



cos(A(i - t'))D(t')dt' . (29) 



Differentiating this equation gives 
8 2 D(t) 



dt 2 



-b 2 D(t) + b 2 A f sm{A(t-t'))D(t')dt', (30) 
Jo 



d 3 D(t) = _ b2 8D(t) 



dt 3 



dt 



b 2 A 2 



= -(b 2 + A 2 ) 



2 .dD(t) 



ot 



cos(A(t - t'))D{t')dt' 



(31) 



As a consequence, the time derivative of D(t) must satisfy 
8D(t) 



dt 



Acos(ft) +Bsm(ft) 



where 



/ = v^ 2 + A 2 



(32) 



(33) 



is the natural oscillation frequency for D(t). From the 
fact that dD(t)/dt vanishes at t = 0, as follows from 
Eq. (|29|), the coefficient A vanishes. At t — the sec- 
ond derivative of D(t) equals — b 2 D(0), so B is given by 
—b 2 D(0)/f. Consequently the derivative of D(t) is given 
by 



^ = -^(0) S in(/*). 



(34) 



Integrating this result and choosing the integration con- 
stant so as to satisfy the initial condition, D(t) is deter- 
mined to be 



D(t) 



A 2 + b 2 cos(ft) 
b 2 + A 2 



D(0) 



(35) 



Expressed in terms of the spin expectation values, note 
that D(0) = -(4)(0), this is 



(5, -/,)(*) 



r& 2 ( l-cos( ft)) 
b 2 + A 2 



(h)(0). (36) 



An analogous calculation of the sum of the spin expecta- 
tion values implies 



(s z + h)(t) 



1 



& 2 (1 - cos(-Ft)) 
6 2 + S 2 



(h)(0), (37) 



where the frequency F = \Jb 2 + S 2 is natural for the 
sum of the spin expectation values. Finally, the sum 
of these results determines the time dependence of the 
S z expectation value, experimentally referred to as the 
"buildup curve" , namely 



(S z )(t) = \b 2 



1 - cos(ft) 1 - cos(Ft) 
b 2 + A 2 b 2 + S 2 



(h)(0)- 



(38) 

Thus it is seen that, for a constant b (specific orientation 
of the vector connecting the S and / spin particles), the 
time dependence of (S z )(t) is governed by a difference 
between oscillating terms. Actually, since lois and u>u 
are both positive (they are the amplitudes of the radio 
frequency fields), the sum frequency E is larger than the 
magnitude of the difference frequency \luis — u)u\. As 
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FIG. 1: Buildup curve for the static Hartmann-Hahn match- 
ing condition u>n = uiis = 60 kHz and a dipolar coupling 
constant bo of 10 kHz. The solid line is for the powder, the 
dashed line for the value of b(8 = n/2) = —10/3 kHz, and the 
dotted line for the value of b(6 = 0) = 20/3 kHz contributing 
to the curve for the powder. 



a consequence, the amplitude b 2 /(b 2 + S 2 ) of the sum 
frequency contribution is smaller than the amplitude of 
the difference frequency contribution, implying that the 
rapid oscillations of the sum frequency are likely lost as 
noise on top of the slower oscillations of the difference 
frequency. _ 

Equation.-d38|) has been obtained by LevittO and 
Zhang et altSTby using a sequence of rotations using ef- 
fective spin operators based on a decomposition of the 
hamiltonian into two commuting parts. 

Experimentally, the matching profile is obtained by 
scanning the rf amplitude of the iS-spin channel and the 
(relatively) longtime average of (S z )(t) measured. This 
should be compared with 



b 2 



6 2 + £ 2 



(/ z )(0). (39) 



This time average is dominated by the difference fre- 
quency term which, if the dominance is sufficient, im- 
plies that the experimental plot of (S z ) t i mc versus u>is is 
Lorcntzian. 



B. Powder average 

In a powder all orientations must be averaged over, so 
the cross polarization time dependence is 



(S z )(t) 



1 



powder 



b 2 (9) 



1 - cos(/(g)t) 
b 2 (0) + A 2 



1 - cos(F(9)t) 
b{9) 2 + E 2 



sm8d6{I z }(0), (40) 



with an analogous expression for the time average. Note 
that since there is no 4> dependence in 6, the average 
over <f> is unnecessary while the frequencies / and F are 
9 dependent through their dependence on b(9). Figure 
1 shows the effects of the angle averaging on the time 



FIG. 2: A comparison between the Hartmann-Hahn matching 
profile for a static CP of a powder sample (solid line, u>n — 
60 kHz and b =10 kHz) and the Lorentzian 2.555/[5.0625 + 
(wis - 60) 2 ] (dashed line). This was fitted at the peak height 
and the width at 1/2 height. 



dependence of (S z )(t), two of the rapidly oscillating con- 
tributions to this powder average are also shown. Com- 
pare the corresponding depolarization curves in Ref. ||. 
The corresponding matching profile, namely a scan of 
the time and powder average (S z ) time powdcr for different 
u>is, is given in Fig. 2. A Lorentzian, picked to have 
the same height and width at half height, is also plot- 
ted. This shows that the powder CP line shape is close 
to Lorentzian, but not quite. 



VI. CROSS POLARIZATION WITH MAGIC 
ANGLE SPINNING 

The orientation of the vector r connecting the posi- 
tions of the two spins is made time dependent by rota- 
tion about an axis n at a frequency of u r . For magic 
angle spinning, the angle between the rotation axis n 
and the direction x of the static magnetic (Zeeman) field 
is arccos(l/v / 3). At that angle the dipolar angle depen- 
dence is split into four time dependent parts, rather than 
the five parts for a general splitting of the order 2 spher- 
ical harmonics. Explicitly this time dependence is given 
by 



b(t) = b* 2 e 



where 



and 



b\e 



V2, 



■tie* 



b 2 e 2 



bi = - — bo cos9 r sm9 r e 1 '^'' 
3 



b 2 = -b sin 2 9 r e 2i ^ 
6 



(41) 
(42) 
(43) 

■At) 



are proportional to the spherical harmonics Y 2 i{9, l 
and Y 22 (9 r ,4> r ), expressed in terms of the angles 9 r , 4> r 
defining the direction between the particles containing 
the spins with respect to the axis of rotation h. A short 
discussion of how this splitting arises is given in Appendix 
B. 
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A. General formulation for magic angle spinning 

The dipolar coupling strength b(t) enters into the inte- 
grodifferential equations (|27j-p!8|) for the cross polarization 
in the form of the product b(t)b(t'). From the explicit 
time dependence given by Eq. (41), the product has a 



total of 16 different terms. Four of these are naturally 
combined into two simple forms involving the time dif- 
ference t — t' , while the remainder have a variety of time 
dependences. In detail, the combination of time depen- 
dences is 



b(t)b(t') = 2\b 1 \ 2 coa(uj r (t-t')) + 2\b 2 \ 2 cos(2uj r (t-t')) 

+b 2 e 2iu r (t+t') + 6 2 e iu,.(t+t') + ( b *j2 e -iu r (t+t') + ( b *}2 e -2i Ur (t+t') 



+hb 2 
+b l b* 2 



,»w r (2t+t') , e iw r (t+2t') 



e iu r (-2t+t') _^ e iw r (t-2t') 



b\b* 2 
+ b\b 2 



e -iui r (2t+t') , e -iui r {t+2t') 
e <w r (2t-t') , e iu r (-t+2t') 



(44) 



In the following it is assumed that only the terms involving the time difference t — t 1 are important, and the 
remainder act to produce a variety of interference effects. Retaining only these two terms, the frequency difference 
formula (27) is approximated by 



dD(t) 



[|6i| 2 cos(w,.(i - t')) + \b 2 \ 2 cos(2w r (t - t'))] cos(A(t - t'))D{t')dt' 
-N- 2 (t) - 2V_i(f) - 2Vi(f) - N 2 (t), 



(45) 



where D(t) is again used as an abbreviation for (S z — 
I z ){t). The integral in this equation is naturally divided 
up into four terms, whose generic definition is 

N n (t) = \b n \ 2 f cos((A + nuo r ){t - t'))D{t')dt' , (46) 



with the obvious definition |6- n | 2 = \b n \ 2 . A solution of 
Eq. d45| ) is required. This can be accomplished as in the 
static CP situation by using successive derivatives of Eq. 
(WEJ , see Sec. But it is necessary to go up to the ninth 
derivative of D(t) before a sufficient number of terms are 
obtained in order to eliminate the N n (tys. This calcula- 
tion is described in Appendix C while an approximation 
is made in the following subsection that provides a sim- 
ple analytical estimation of the cross polarization with 
magic angle spinning. 



B. Independent Sideband Approximation 

The idea here is that, near a Hartmann-Hahn matching 
condition, only one of the -/V„(t)'s is of importance, that 
is, of appreciable magnitude. As a consequence, D(t) 
can be considered as a sum of four terms, each of which 
is determined by a separate equation. With a similar 
approximation for the T, based G(t) = (S z + I z )(t), the 
total cross polarization is described as a sum of eight 
terms. 

To implement this idea, D(t) is written as the sum 
D{t) = D_ 2 (i) + + Di(i) + D 2 {t), (47) 



each term of which is to be approximated as the solution 
of 



= -16, 



| 2 / cos((A +nu r )(t-t'))D n (t')dt'. 
Jo 



at 

(48) 

A decoupling of Eq. (|45|) has thus been accomplished. 
Since the equation for D n (t) has the same structure as 
that of Eq. (|29|), so its solution is, according to Eq. (p5q), 



D n (t) 



(A + nuj r ) 2 + \b n \ 2 cos(f n t) 



\b n \ 2 + (A + nuj r ) 2 
with an effective frequency given by 

/„= VK\ 2 + (A + nu; r y 



D(0), 



(49) 



(50) 



On the basis that only one D n (t) has a significant con- 
tribution when near a matching condition, it has been 
assumed that D n (0) is well estimated by -D(O). Adding 
in the corresponding approximation for (S z + I z ){t) and 
doing some minor algebra, the complete time dependence 
of (S z ){t) is estimated by the sum of eight terms 



1 



E 



K\ 



n=— 2,-1,1,2 



1 — cos(/„i) 
|&„| 2 + (A + nuv) 5 



1 — COS^nt) 



W(0). (51) 



|6„| 2 + (£ + n^.) 2 
Here the frequencies for the sum terms are given by 
F n = ^/\b n \ 2 + (E + nto r ) 2 . 



(52) 
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A powder average of this result is easily calculated since 
the |6„| 2 's depend only on the single angle 9 r . The time 
average of Eq. ( |5l| ) is obtained by merely dropping the 
cosine terms, to give an estimation of the CPMAS match- 
ing profile. 

b =10kHz 

0.5-1 b„ = 1 kHz 



0.4- 




30 40 50 60 70 80 90 



co ls [kHz] 

FIG. 3: The CPMAS matching profile for a powder sample 
with u) r = 10 kHz and Un = 60 kHz. The solid line is for a 
dipolar coupling constant 6o of 10 kHz and the dashed curve 
for a dipolar coupling constant of 1 kHz. 

Figures 3 and 4 illustrate the typical matching pro- 
files that can be attained. Figure 3 exemplifies what 
is standardly observed in an experiment, namely four 
peaks at the shifted Hartmann-Hahn matching condi- 
tions u>is = + nuj r for rt=-2,-l,l and 2. These are 
due to the set of first terms in the sum of Eq. ( p3l| ) 
[with the cosine term dropped to reflect the time aver- 
age], associated with the conditions A + nuj r = 0. In 
such a situation the second set of terms is small because 
T,-\-nuj r A + nLu r . But if 2cu r > ujh then lou — 2u> r < 
and the matching condition A + 2uj r requires a negative 
uj\s- This is not experimentally possible (wis measures 
the strength of an rf field) so this matching condition 
cannot be attained and the peak in the matching profile 
is absent. However, in this case, one of the second set 
of terms in Eq. (|5l| ) can have a small denominator, so a 
peak in the matching profile arises from that term. This 
gives a negative (S z )(t) so the peak is negative. Figure 
4 illustrates this case. The structure is emphasized by 
choosing a small uu while a large uj r is needed to allow 
the matching condition, u>is = 2ui r — u>u to be attained 
for positive u>is- MeierQ has a figure of experimental data 
illustrating this kind of behavior, but for the 3-spin sys- 
tem CH2. An even more interesting profile, which does 
not seem to have yet been observed experimentally, is 
when lo\i is even smaller, so that two matching condi- 
tions are from the second term, namely £ — 2uj r = and 
£ — uj r = 0, while the other two matching conditions are 
A — Lj r = and A — 2cu r . The result is a matching profile 
with two negative and two positive peaks. This is illus- 
trated by the second profile shown in Fig. 4. It is clear 
that there can be at most only two negative peaks in a 




10 20 30 40 



co ls [kHz] 

FIG. 4: The CPMAS matching profiles for a powder sample 
with a dipolar coupling constant bo of 10 kHz. The solid line 
is for urn = 12 kHz and ui r = 10 kHz while the dashed line is 
for uj\i = 8 kHz and uj r — 11.34 kHz. 



matching profile. 



VII. DISCUSSION 

The special properties of spin-1/2 systems and the spe- 
cial nature of the hamiltonian, Eq. (ra), have been used 
to give a formally exact method of calculating the cross 
polarization buildup curves and matching profiles for an 
isolated pair of spins. For the static case (no sample spin- 
ning), this can be expressed in closed form, Eq. (pq), if 
the interspin direction is fixed, while a powder average 
can be easily calculated numerically, Eq. ■— fio|) . These 
results have also been obtained by Levittliil and Zhang 
et a£3 using a sequence of rotations based on a partic- 
ular division of the hamiltonian into commuting parts. 
It is of note that the time dependence of the cross po- 
larization is purely oscillatory if a single orientation of 
the dipolar coupling strength is involved. It is the sen- 
sitive dependence of the frequency of this oscillation on 
the orientation, which causes an enormous destructive in- 
terference in a powder sample. The decay-like nature of 
the envelope of the oscillations in a powder sample, see 
Fig. 1, arises from this destructive interference rather 
than from any spin relaxation (or spin diffusion), which 
of course is not included in the hamiltonian dynamics of 
an isolated system. 

Magic angle spinning is formally easy to include, it is 
just a case of recognizing the correct time dependence of 
the dipolar coupling strengths in the integrodifferential 
equations ( p7| ) and ( p8| ) . But how to obtain a general so- 
lution of these equations is not clear. Two simplifying ap- 
proximations have been made. First, only those terms in 
the product of the dipolar coupling strengths, Eq. (jij), 
are kept that depend on the time difference. It is assumed 
that the other time dependences lead to minor interfer- 
ence effects but that has not been explored. A method 
of solving the equations with the above single approxi- 
mation is given in Appendix C, and been explored with 
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some calculations. The results found, which are not pre- 
sented here, indicate that the solutions are qualitatively 
in agree ment with those using the further approximation 
of Sec. VI B . This second approximation assumes that 
only one term in Eq. ( fl5| ) is of importance in a range 
of wis values. An analytic expression, Eq. (|5l]), for the 
cross polarization is then obtained for fixed interspin ori- 
entation, and it is in a form in which a powder average 
can easily be carried out numerically. Its form is very 
simple and it is immediately seen how the cross polariza- 
tion is significant when the standard matching conditions 
A + mjj r = are met. It is also seen that there is another 
set of matching conditions, E + nuj r = which can lead 
to an enhancement of the cross polarization, but now 
with (S z ) and (I z }(0) having opposite signs. There can 
be at most two negative peaks in the matching profile, 
but there are always four physically possible matching 
conditions for a given u>n and u) r . 

The method presented here, with the approximations, 
has the advantage that the qualitative behavior of ex- 
perimental buildup curves and matching profiles can be 
given simple explanations. The quantitative connection 
between the dipolar coupling strength and the size and 
width of the matching profiles, and the oscillatory nature 
of the buildup curves is formally available. How well that 
can be done experimentally can also depend on the real 
physical system. The spin pair is never truly isolated, 
so small interactions can cause some more interference 
effects, and the surroundings afford relaxation. Such ef- 
fects can surely influence the shape of the buildup curves, 
contributing to a more rapid closure of the envelope of 
the oscillations as time progresses. It would also affect 
the height of the matching profile peaks, but hopefully 
not perturb too much the position of the peaks. 

It is interesting that the thermodynamic method can 
give the same formulas for the A-based peaks in the 
matching profile as the independent sideband approxima- 
tion, provided an aj r -dependent transformation is used 
to get a time independent rotating frame hamiltonian. 
However, this must be done one peak at a time whereas 
the present method gives all four peaks predicted by the 
same formalism and there is no question about any fur- 
ther frame transformation. 

Of the previously introduced methods that we are 
aware of £er interpreting cross polarization, the one used 
by LevittO and Zhang et afca seems the closest to being 
fundamentally based. But it does not seem to have been 
extended to include the magic angle spinning case. It 
is also not clear to us how that might be incorporated 
into their formalism. We see the formula introduced by 
MiillerO, and extensively used for fitting experimental 
buildup curves, as a very useful fitting tool. But since 
it inherently has exponential decays built into it, care 
should be taken when interpreting the parameters found 
in terms of physical quantities, since no relaxation is fun- 
damental to the explanation of the buildup curves. We 
hasten to add that, as we indicated above, in a real sys- 
tem the spins are not truly isolated so relaxation effects 



are present and these can influence the detailed shapes 
of the experimental curves, but not their fundamental 
behavior. 
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APPENDIX A: EIGENVALUES AND 
EIGENVECTORS OF H 

This lists the eigenvectors and eigenvalues of the hamil- 
tonian defined in Eq. (^) for time independent b and for 
fixed orientation. The basis \mn), m,n — ±1/2 used for 
this calculation is that of the simultaneous eigenvectors 
of S z and I z , namely 

S z \mn) = m\mn); I z \mn) = n\mn). (Al) 

As a notation similar to Zhang'sEl, the quantities 

, X 



A = y/b 2 + A 2 ; 



VP 



(A2) 



are introduced. Then the eigenvectors, labelled with 
their respective energy eigenvalues are: 



l^(l/2)w a ) = 
IV>-(1/2)u>a) = 
1^(1/2)^) = 



b\l/2, -1/2) + {uj a - A)| - 1/2, 1/2) 



(l/2)u 





- {uo A - A) 2 


b\- 1/2,1/2)- 


(c A -A)|l/2,-l/2) 




- {lo a - A) 2 


ft|l/2,l/2) + (w 


s - E)| - 1/2,-1/2) 




- (w s - E) 2 


61-1/2,-1/2) 


-( W s -E)|l/2,l/2) 




- (w s - E) 2 



It is be noted that these are NOT the same as the eigen- 
values and eigenvectors of Zhang's H A and/or -ff s , aris- 
ing from his decomposition of the hamiltonian of Eq. (|^), 
even taking into account the difference in the labelling of 
the coordinate axes. 

The limits of these eigenvectors when 6 — > and A — > 
or E — * give respectively the eigenvectors of the rf 
hamiltonian and of the truncated dipolar hamiltonian. 
In particular, as b — > with A > 0, the expansion of uj a 
is 



to 



A 



A A 1 



A 2 



A- 



2A 



(A3) 
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while u> A — > b as A — > 0. As a consequence, the first listed eigenvector has the limits 
I 

Ml/2, -1/2) asfo^0forA>0 
|Vi/ao,A) < I - 1/2, 1/2) as fo^ for A < (A4) 

I (1/V2) [|l/2,-l/2) + | - 1/2, 1/2)] as A 0. 

The limits of the other eigenvectors are obtained in an analogous manner. 

APPENDIX B: MAGIC ANGLE SPINNING 

A vector based derivation is given for the effect of spinning the orientation dependent dipolar coupling strength. 
The rotation of any unit vector r about an axis n by an angle x produces the rotated vector 

f = hh-f + (f — hh-f) cosx + nxf sin^- (Bl) 

The component of the rotated vector along the Zeeman field direction x is then 

cos 9 = x-r = x-hh-f + (x-f — x-hh-f) cosx + xXn-fsmx- (B2) 

If the unit vector f is expressed in a right handed coordinate frame based on the rotation axis h and the x — h plane, 

x — hh-x . hxi; 

r = n cos 9 r H z — sin 9 r cos <fi r H - — sin 9 r sin </> r , (B3) 

y/\ — (h-x) 2 y/l — (h-x) 2 

then 

cos 9 = x-hcos9 r + \J\ — (h-x) 2 sin9 r cos(</v + x)- (B4) 

The angle dependence of the rotated dipolar coupling strength is then expressed by 

cos 2 9 -1/3 = (h-x) 2 cos 2 9 r - 1/3+ [1 - (h-x) 2 ]sm 2 9 r cos 2 (<j) r + X ) 

+2n-x\Jl — (h-x) 2 cos9 r sin# r cos(0 r + x). (B5) 

The trigonometric identity 

cos 2 ((j )r +x) = i[cos2(0 r + X ) + l] (B6) 

simplifies this expression. Finally, if h-x is chosen as l/y/3 and the angle x is se t equal to the w r t, the time dependence 
of the dipolar coupling strength is 

/ ^ \ ^ 2^/2 

bo I cos 2 9 — - ) = 6o- sin 2 fl r cos[2(0 r + uj r t)] + bo » cos 9 r sinfl r cos(0 r + cj r t) 

V tj J o o 

= b 2 e T +O x e + Die + o 2 e r , (B7J 

where &i and &2 are given explicitly in Eqs. ( [l2| ) and (Jl3|). 

It is noted that the sum of the squares of the magnitudes of these parts of the coupling strengths, when averaged 
over all angles, preserves the total square of the coupling strength, namely 

^(cos 2 ^- i) 2 sin 9d9d<p = J J 2\b x \ 2 + 2\b 2 \ 2 sin 9 r d9 r d<j) r 

= ^ = 1^6 2 + ^U 2 . (B8) 
45 135 135 v ; 

Thus the contributions of the b\ terms are twice those of th e 62 terms. 

APPENDIX C: COMPLETE SIDEBAND the second derivative of N n (t) is 



CALCULATION 



d 2 N n (t) 2 dD(t) 



dt 2 dt 
This appendix gives a method for carrying out the ex- 
act calculation of D(t) from Eq. pa). For this purpose, 



= K\ ~ ( A + ntu n ) 2 N n (t). (CI) 
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Then the third derivative of D(t) can be expressed in 
terms of the first derivative and the iV„(i)'s, that is 



d 3 D(t) 
dt 3 



Y J K\ 2 ^^+Y.^+™n?N n {t). (C2) 



Continuing in this way, the higher odd-ordered deriva- 
tives of D{t) can be calculated in terms of lower ordered 
derivatives and the iV„(£)'s. The first 4 equations ob- 
tained in this way (up to the 7 th order derivative) can 
be used to solve for the N n (t)'s in terms of the deriva- 
tives of D{t). Substituting these into the equation for the 
9 th order derivative produces a differential equation for 
D(t) of 9 th order, but with constant coefficients. Such an 
equation has 9 solutions of the form D(t) = e xt , whose 
proper combination gives the time dependence of D(t). 

To simplify this calculation, it is convenient to intro- 
duce the notation D^ m \t) for the mth derivative of D(t), 
and the sums 



SS m = ^2(A + mj r ) m \b n 



(C3) 



Then the needed equations for the higher odd-ordered 
derivatives of D(t) are 

D i5 \t) = -SSoDW +SS 2 dW -^2(^ + nto) 4 N n (t) 

n 

D (7 \t) = -SS a D (5) + SS 2 D {3) - SS 4 D {1) 
+ ]T(A + r^) 6 A n (i) 

■n 

D {9) (t) = -SS D {7) + SS 2 D {5) - SSiD {3) 

+55 6 D (1) -^(A + mj) 8 7V„(t). (C4) 

n 

A formal inversion of the matrix of frequency powers 
multiplying the iV n (£)'s from the 4 odd-ordered deriva- 
tives £>W(i) to DW(t) allows the iV n (t)'s to be expressed 
in terms of the £>( m )(i)'s. Substitution into the D^(t) 
equation gives the desired differential equation for D(t). 
This involves only odd-ordered derivatives of D(t). If e At 
is substituted in for D(t), this gives a fourth order poly- 
nomial in A 2 times a factor of A. The polynomial has 4 
roots which must be negative since the spin dynamics is 
governed by the quantum Liouville equation (and thus 
the spin system must oscillate without any relaxation), 
and there is also the zero root from the multiplicative 
factor of A. The 4 eigenfrequencies fi m (square roots of 
minus the polynomial roots) and the zero frequency de- 
termine the time dependence of D(t). A further simplifi- 
cation arises from noticing that N n (0) = 0, which implies 
that all odd derivatives of D(t) vanish at t = 0. Thus 
D(t) must be an even function of t and can be expressed 



as the sum 

4 

D(t) = A + A ™ cos(O m t). (C5) 

m— 1 

The expansion coefficients A m (m = • • • 4) still need to 
be determined. 

The coefficients A m are determined from the even- 
ordered derivatives D^(t) evaluated at t = 0. These 
can be found directly by differentiating Eqs. (^5|), ( |C^ ) 
and (p4p, and evaluating them at t — 0. Specifically, the 



5 required relations 


are 




A a +^A m = 

m 


D(0) 




— ^ ^m A m = 

m 




= -SS Q D(0) 


m 




= (SS 2 + SS 2 )D(0) 


— ^ £l & m A m = 


D {6 \0) 


= -[SS (SS% + 2SS 2 ) 


m 


-i 


■SSa]D(Q) 


^ ^m A rn = 




= [SS%(SS% + 3SS 2 ) 



+2SS SS 4 + SS% + SS e ]D(0). (C6) 

These equations are to be solved for the expansion coef- 
ficients. 

This determines the difference (S z — I z )(t) = D(t) of 
the expectation values of the two spin operators. A cal- 
culation of the sum (S z +I z )(t) = G{t) involves the same 
treatment, but with the sum frequency £ replacing the 
difference frequency A in all the above formulas. The 
cross polarization is given by one half the sum of these 
two quantities 



{S z ){t)= l -[D(t)+G{t)]. 



(C7) 



Note that G(0) = -D(0) = (I z )(0). To calculate the 
powder average CPMAS, the above computations need 
to be carried out for every possible orientation of the 
vector pointing from one spin to the other. A bonus of the 
approximation of retaining only the terms |&i| 2 and \b 2 \ 2 
in the product of the time dependent coupling strengths, 
Eq. (Q) , is that these quantities are independent of the 
angle <p r in the description of this orientation, so only the 
9 r average needs to be done. Thus the powder CPMAS 
is given by 



(S 



z I powder 



it) 



[D(t) +G(t)]d cos ( 



(C8) 
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